Interstitial Mn in (Ga,Mn)As: Binding energy and exchange coupling 
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We present ab initio calculations of total energies of Mn atoms in various interstitial positions. 
The calculations are performed by the full-potential linearized plane-wave method. The minimum 
energy is found for tetrahedral T(As4) position, but the energy of the T(Ga4) site differs by only 
a few meV. The T(Ga4) position becomes preferable in the p-type materials. In samples with 
one substitutional and one interstitial Mn the Mn atoms tend to form close pair with antiparallel 
magnetic moments. 

We also use the spin-splitting of the valence band to estimate the exchange coupling Jpd for various 
positions of Mn. The exchange parameter is the same for the substitutional and for the T(As4) 
position and it is somewhat smaller in the case of the T(Ga4) position. 

PACS numbers: 71.15.Mb, 71.15.Nc, 75.50.Pp 
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I. INTRODUCTION 

Diluted magnetic semiconductors, such as (Ga,Mn)As, 
are important materials combining ferromagnetic behav- 
ior with a sensitivity to doping characteristic for semicon- 
ductors. This leads to their interesting physical proper- 
ties and makespossiblc their applications e.g. in spin- 
electronics H H. Although (Ga,Mn)As has bee n ex- 
tensively studied in the last years, some aspects of the 
incorporation of Mn into the crystal lattice still remain 
unclear. It was assumed that in the well defined samples 
Mn simply substitutes for the host cation. Only recently 
it was pointed out 0, 01 that the presence of interstitial 
Mn may explain some peculiar properties of (Ga,Mn)As 
such as the low doping efficiency of Mn acceptors p| . The 
main reason for considering the Mn atoms on the inter- 
stitial positions was that they act as double donors and 
partly compensate the Mn acceptors in the substitutional 
positions. 

At the same time, channelling Rutherford backscat- 
tering experiments proved that a large fraction of Mn 
atoms indeed occupies the interstitial positions (Mnint) 
in the as-grown samples 0. A correlation between the 
removal of Mn interstitials and increase of the conductiv- 
ity, the Curie temperature, and saturation magnetization 
has been found. 

In addition to the self-compensation effect, the inter- 
stitial Mn atom also reduce the number of local moments 
that participate in the ferromagnetic state. This was ex- 
plained by pairing of Mnjnt with the Mn atoms in the 
substitutional positions (Mnca) due to their Coulomb 
attraction At the bonding distance, the antiferro- 
magnetic superexchange within the pair is assumed to 
outweigh the hole mediated ferromagnetic exchange. As 
a result, the moments of the paired Mn atoms have op- 
posite directions and the pair as a whole has no magnetic 
moment. 

The ab initio studies of the interstitial Mn showed the 
differences in the electronic structure of the Mn atom in 
the interstitial and substitutional positions 0,01 and also 
a possible reaction path for the incorporation of Mn into 



the GaAs lattice [g- It was also found that the increase 
of the lattice constant of Gai_xMnxAs with increasing 
content of Mn is partly due to the presence of Mn in the 
interstitial positions 9]. 

Recently, Blinowski and Kacman [Tol investigated 
the coupling of the local moment on the interstitial Mn 
to the spin of the carriers in the valence band. Using a 
simplified tight-binding model they showed that the local 
moment on Mn atom in the tetrahedral position with four 
Ga neighbors, T(Ga4), is effectively decoupled from the 
spin of the holes in the valence band. 

They claimed that the coupling is weak because, in 
addition to the reduced number of the holes, also the 
coupling constant Jpd is small for the interstitial Mn in 
the T(Ga4) position. As a result, the local moment is not 
subject to the ferromagnetic coupling with the moments 
around it. This opens the way for the AF exchange to 
be irnportant in the Muint-MnQa pair, as anticipated in 
Ref. II. 

There are, however, several open questions concerning 
the magnetic interactions of Mn in the interstitial posi- 
tions. First of all, the spin-polarized band structures 0| 
did not show any indication of different values of Jpd for 
substitutional and interstitial Mn, at least for the T(As4) 
positions. In addition, the Mnint-MuQa pair as a whole is 
expected 0,113 to have only a small magnetic moment, 
but the exchange coupling with MuQa ^ being uncom- 
pensated by the contribution of Mnint ^ should strongly 
polarize the holes. 

That is why we performed a more detailed study of 
the interstitial Mn and its spin interactions. We use the 
density-functional, full-potential linearized augmented 
plane wave calculations (FPLAPW to obtain the 

electronic structure of (Ga,Mn)As with Mn atoms in var- 
ious crystallographic positions. The calculated total en- 
ergies are used to compare different positions of the in- 
terstitial Mn and to estimate the strength of the MuQa- 
Muint pair interactions. The splitting of the valence band 
for the majority- and minority-spin electrons is used to 
compare the corresponding values of the exchange pa- 
rameters Jpd- 
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II. Mn IN VARIOUS INTERSTITIAL 
POSITIONS 

We compare the total energies for three interstitial po- 
sitions of Mn in GaAs. There are two inequivalent tetra- 
hedral positions in the zinc-blende structure of GaAs, 
T(As4) and T(Ga4). They are surrounded by four As 
and Ga atoms, respectively. The (unrelaxed) distances 
of these nearest neighbors are equal to the length di of the 
Ga-As bond, i.e. to the Mn-As distance for the substitu- 
tional MuQa. The local arrangement around the intersti- 
tial Mn is, in contrast to MuQaj characterized by another 
six close neighbors at the distance ^2 = \/9/8 • di « 
1.155(ii. In the hexagonal position, the interstitial Mn 
has three Ga an d three As atoms at the same distance 
dhex — 1/11/12 • di « 0.957c?i, and no other close neigh- 
bors. 

We represent GaAs with Mujnt by hexagonal 
Gai2Asi2Mn supercells. The c-axis coincides with the 
body diagonal of the conventional cubic unit cell. The 
symmetry of the cell does not change if we shift the Mn 
interstitial along the c-axis from T(As4) to hexagonal 
and T(Ga4) positions. This makes possible to perform 
all calculations under the same conditions. 



Table I: Total energy, Etot, and total spin, Stot, of the unit 
cell of Gai2Asi2Mn with various interstitial positions of Mn. 
Ga9Zn3Asi2Mn samples are used to simulate p-type materi- 
als. 



Sample 


Mn position 


Etot (eV) 


Stot 


Gai2Asi2Mn 


T(As4) 


ground state 


1.665 


Gai2Asi2Mn 


T(Ga4) 


+ 0.005 


1.555 


Gai2Asi2Mn 


hex. 


+ 0.522 


1.519 


Ga9Zn3Asi2Mn 


T(As4) 


+ 0.063 


2.331 


Ga9Zn3Asi2Mn 


T(Ga4) 


ground state 


2.176 


Ga9Zn3Asi2Mn 


hex. 


+ 1.320 


2.136 



The results are summarized in Table I. Surprisingly, 
the binding energy of the interstitial Mn does not depend 
much on its nearest neighbors. The interstitial Mn has 
minimum energy in the T(As4) position. However, the 
difference of the total energies obtained for Mn in T(Ga4) 
and T(As4) positions is of order of a few meV and can be 
neglected in practice. This means that, without interven- 
tion of other charged defects, the Mn interstitials can be 
found with an almost equal probability in either T(As4) 
or T(Ga4) position. The total energy corresponding to 
the hexagonal interstitial position of Mn is remarkably 
higher and represents a barrier w 0.5 eV separating the 
tetrahedral positions. Calculations with similar results 
have been recently done by Boguslawski 13]. 

The weak influence of the nearest neighbors on the 
interstitial Mn in T(As4) and T(Ga4) positions can be 
observed also in the densities of states in Fig. 1. The 
total DOS and the distribution of the Mn d-states are al- 
most identical. Also the DOS of Gai2Asi2Mn with Mn in 
the hexagonal position is quite similar and differs mainly 
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Figure 1: Spin-polarized densities of states for Gai2Asi2Mn 
with an interstitial Mn atom in T(As4) (upper panel), T(Ga4) 
(middle panel), and hexagonal position (lower panel). The 
contribution of Mn d-states is indicated by hatched area. 



by the overlap of the valence band for the majority-spin 
electrons (|) with the conduction band for the minority- 
spin (J.) electrons. This tendency to close the gap in the 
electron spectrum correlates well with the increase of the 
total energy. 

The lattice relaxation around the Mn impurity is not 
very important for the substitutional Mn as shown by 
Mirbt et al. 0|. This is not the case of the intersti- 
tial Mn which tends to expand the crystal lattice j^. 
The addition of Mn into either T(As4) or T(Ga4) posi- 
tions results in a remarkable repulsion of the nearest and 
next-nearest neighbors. The positions of the more distant 
atoms in the supercell are changed much less. The energy 
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gain due to the relaxation is in both cases approximately 
20 meV and it does not change the ordering of the total 
energies for the T(As4) or T(Ga4) positions. The lattice 
relaxation around Mn in the hexagonal position has not 
been considered because of its minor importance. 

In the case of the T(As4) position, the distance of the 
four nearest As neighbors was found to increase by « 
1.5 % from 2.45 A to 2.47 A. The distance of the six 
next-nearest Ga neighbors increases by « 0.7 % from 
2.835 A to 2.88 A. 

The relaxation around Mn in the T(Ga4) position is 
larger. The nearest Ga neighbors are pushed to the 
distance 2.515 A and the relaxed distance of the next- 
nearest neighbors increases by « 0.5 %. The enhanced 
relaxation of the nearest neighbors in this case is due to 
the Coulomb repulsion between Mn and Ga atoms which 
are both positively charged. 

In all cases, there are two electrons in the conduction 
band, i.e., the Mnint always acts as a double donor. The 
electrons in the conduction band are almost completely 
spin-polarized. They accumulate in the minority-spin 
conduction band, so that the total spin of the cell (i.e., 
the spin per Mn) is reduced to ~ 3/2 in accordance with 
our previous calculations 

For comparison, we performed the same set of calcula- 
tions also for hypothetical GagZn3Asi2Mn crystals with 
Zn atoms substituted at the sites most distant from the 
interstitial Mn. The presence of Zn has only a little effect 
on the density of states, but the material is converted into 
the p-type with one hole in the valence band. In this case, 
the T(Ga4) turns to be the stable interstitial position 
of Mn. The energy difference between the T(As4) and 
T(Ga4) positions, approximately 60 meV, is high enough 
for a preferential occupation of the T(Ga4) position at 
typical growth conditions. The remarkable increase of 
the energy of the hexagonal position indicates that the 
mobility of the interstitial Mn may depend on the type 
and degree of the doping. 



III. BINDING ENERGY AND EXCHANGE 
COUPLING OF A Mn - Mn PAIR 

The hexagonal unit cell used in Sect. II is well suited 
also for the study of the Muint-Mnca pair interactions. 
We consider a hypothetical GaiiMnAsi2Mn crystal with 
one substitutional and one interstitial Mn in the unit 
cell. We consider three positions of the interstitial Mn 
shown in Fig. 2. The positions T(Ga3Mn) and T(As4), 
denoted (a) and (b) in Fig. 2, are located on the c-axis 
of our unit cell and correspond to the initial stage of 
dissociation of the Mnca-Muint pair. On the other hand, 
the T(As4) position denoted (c) is very close to MuQa 
and is - together with T(Ga3Mn) - a candidate for the 
ground state of Muinf 

In the T(Ga3Mn) position, Mujut and MuQa are the 
nearest neighbors at the distance di. We performed 
the density-functional calculations with two initial con- 



ditions, corresponding to parallel and antiparallel align- 
ment of their local moments. In both cases, the self- 
consistent procedure converges to a locally stable solu- 
tion without changing the initial alignment of the local 
moments. The resulting total energies are given in Ta- 
ble II. The antiparallel alignment is energetically more 
favorable than the parallel alignment, in a good corre- 
spondence with the expectations ^ .11] . The coupling is 
strong enough so that the AF state of the Muint-Mnca 
pair is stable with respect to the thermal fluctuations. 

The Table II shows also the local spins on Mnca and 
Muint atoms defined as integrals of the spin density 
over the corresponding atomic spheres. Although these 
quantities are not directly related to the size of the 
observable local moments, we can see that the local 
moments of MuQa and Muint are comparable and that 
the total magnetic moment of the pair in the ground 
state is much smaller than the magnetic moment of a 
single Mn. 

Assuming that the exchange coupling between more 
distant Mn atoms from different unit cells is much smaller 
than the exchange interaction J{di) within the closest 
pair we can estimate its strength from the difference 
AE{di) = E^^{di) - E^iidi) « ~2J{di)S^ of the to- 
tal energies. Using saturated values S = 5/2 for both 
local moments we obtain the lower estimate for J(di), 
namely J{di) w -26 meV. 

The coupling between Mni„t and MuQa remains anti- 
ferromagnetic also if Muint moves to any of the adjacent 
T(As4) positions (b) and (c) in Fig. 2. The energy dif- 
ference between the parallel and antiparallel alignment of 
the local moments decreases with the increasing distance 
of the Mn atoms, as expected for the superexchange. It 
remains almost the same for the nearest and close next- 




Figure 2: Configurations of the Mnoa-Mnint pair in the (110) 
plane of GaAs crystal [irI |: (a) Mnint in T(Ga3Mn) position 
is the nearest neighbor of Mnoa; (b) Mnint in T(As4) position 
representing a partially dissociated pair at a doubled distance; 
(c) Mnint in the T(As4) position closest to Mnoa. 
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Table II: Total energy of GaiiMnAsi2Mn unit cell and the 
spin assigned to Mnint and Mnoa atoms for various configu- 
rations of the Muint-MuGa pair shown in Fig 2. The config- 
uration (c) is treated separately because of its different sym- 
metry. 



Pair 


MuGa - Mnint 


Total energy 


Spin 


Spin 


config. 


distance (A) 


(eV) 


of Mnoa 


of MUint 


(a) TT 


2.443 


+ 0.324 


1.678 


1.617 


(a) n 


2.443 


ground state 


1.778 


-1.531 


(b) TT 


4.886 


+ 0.502 


1.930 


1.616 


(b) Ti 


4.886 


+ 0.330 


1.899 


-1.558 


(c) TT 


2.835 


+ 0.303 


1.774 


1.657 


(c) Ti 


2.835 


ground state 


1.842 


-1.549 



nearest pairs and it is reduced approximately to one half 
of AE{di) for the doubled distance corresponding to the 
configuration (b). 

It should be pointed out, however, that the density- 
functional calculations tend to overestimate the strength 
of the exchange coupling. The reason for this is that the 
exchange splitting of the Mn d-states, i.e. the separa- 
tion of occupied majority-spin and empty minority-spin 
states on the energy scale levels is systematically under- 
estimated. According to our calculations, the exchange 
splitting Sd{l) — EdCT) deduced from the spin-polarized 
spectral distribution of Mn d-states ranges from 2 eV to 
3 eV. This is roughly one half of the realistic estimate for 
the exchange splitting in Mn (cf. e.g. Ref. ^3)- Cor- 
respondingly, the above given value of J{di) should be 
divided by four. In this way, we end with a value close 
to the result of Ref. [ill- 

Finally, we estimate the binding energy of the Muint- 
MuQa pair. We compare the energies corresponding to 
the (a) and (b) configurations from Fig. 2, both in the 
magnetic ground state with the antiparallel alignment of 
the local moments. The partially dissociated pair has a 
higher energy and the energy difference i(d2) — E-^ i(di) 
is approximately 0.33 eV. It is, however, only the lower 
estimate for the binding energy of the Mnjnt-MnGa pa-ir 
because the dissociation of the pair is far from being 
complete in our periodic model. Nevertheless, even the 
value obtained here indicates that the Mn interstitials are 
strongly attracted by the Mnca atoms. As long as the 
concentration of Mni„t is lower than the concentration 
of MuQa we can assume that most of Mnjnt atoms are 
involved in the pairs and that the blocking mechanism 
proposed in Ref. works. 



IV. EXCHANGE INTERACTION OF 
INTERSTITIAL Mn AND HOLES 

Due to the hybridization with the spin-polarized 
Mn d-states, also the distribution of the states derived 
from the GaAs valence band depends on the spin. In 
particular, the valence-band states for the majority-spin 



Table III: Spin splitting AE^, — Ev{]') — Ev{l) for ferromag- 
netic state of (Ga,Mn)As. 



Sample 


Mn positions 


A 7^ 


7 

'J pa 






[e\ ) 


[^G V nm ) 


LTai2MnASi2 


i (AS4j 


U.DZ 1 


n 1 /I 

U. 14 


ijai2iVlnAsi2 


i (Lja4 ) 


u.oio 


nil 
U.li 


Gai2MnAsi2 


hex. 


0.653 


0.14 


GaieMnAsie 


T(As4) 


0.502 


0.14 


GaieMnAsie 


T(Ga4) 


0.429 


0.12 


Gai4Mn2Asi6 


2 X Mnca 


1.021 


0.15 


Gai4Mn3Asi6 


2 X Mnca + T(As4) 


1.328 


0.13 


Ga3oMn3As32 


2 X Mnca + T(As4) 


0.743 


0.14 



electrons hybridize with the occupied d-states and are 
pushed to higher energies. The minority-spin states 
are, on the other hand, pushed down due to their 
hybridization with unfilled d-states. 

This effect is formally described by the Kondo ex- 
change interaction between the local spins at sites Ri 
occupied by Mn and the spin density s(r) due to the 
itinerant holes p^ . 

Hint — Jpd^'^^i ■ ^iRi)- (1) 

i 

The exchange parameter Jp^ characterizes the strength 
of the coupling. Within the mean-field theory, the Kondo 
exchange interaction results in the splitting of the valence 
band edge Ey. The splitting AEy is proportional to the 
size S of the local spins and to the concentration x of 
magnetic ions, 

4-x 

AE, = E,{1) ~ ^,(t) = -^JpdS (2) 

assuming spin ^ for the holes, a is the lattice constant. 
We use Eq. (|2Jl to determine the exchange parameter Jpd 
from our spin-polarized band structures. The results of 
the calculations for the unrelaxed geometries presented 
here and in Ref. are summarized in Table III. 

The resulting values of Jpd are overestimated from the 
same reasons as discussed above in Sect. III. Being di- 
vided by a factor of two (the reduction factor for the en- 
ergy nominator) , they approach the realistic values |0| . 
In this work, however, we concentrate on the comparison 
of Jpd for different geometries. The overestimate of Jpd 
due to the reduced band gap is expected to be similar in 
all cases and not very important in this respect. 

In Fig. 3, we plot AEy against x to visualize the lin- 
ear increase of the band splitting with the concentration 
of Mn. The most noticeable result is that all points in 
Fig. 3, obtained for substitutional Mn as well for Mn in 
the interstitial T(As4) and T(Ga4) positions lie around 
the same line. This means that the value of Jpd does 
not depend much on the position of the Mn atom in the 
lattice. 
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simplified tight-binding models mentioned in the Intro- 
duction. 

In addition, the effective 'annihilation' of Mnca due to 
pairing with Mnint is not restricted to the closest pairs 
with Muint in the T(Ga3Mn) position, but it works as 
long as the Mn atoms are close one to another and the 
exchange coupling in the pair remains antiferromagnetic. 
V. CONCLUSIONS 



0.0 I ' ' 1 

0.00 0.05 0.10 0.15 

Concentration of Mn 

Figure 3: Spin splitting of the valence band as a function of 
Mn concentration. 



This result is not very surprising in the case of the 
substitutional Mn and T(As4) position. In either case, 
Mn is surrounded by four As atoms at the same distance 
di and the same degree of the hybridization of the valence 
band states with Mn d-states can be expected. 

The results for the T(Ga4) position are not simple to 
interpret. The value of Jpd for both Gai2MnAsi2 and 
GaigMnAsie is smaller than but comparable with the 
result for the substitutional Mnca site. Assuming for 
simplicity that only the hybridization of the Mn d-states 
with p-statcs of the nearest neighbors is relevant for the 
spin splitting AEy, the values of Jpd simply reflect the 
orbital composition of the valence band. It is well known 
that the top of the valence band in GaAs, as well as 
in other III-V semiconductors, is composed of both an- 
ion and cation p-states, in the proportion approximately 
c^" 3/4, cf'' X 1/4 113. This means that the same 
proportion (roughly 3:1) should be expected also for the 
exchange parameters Jpd parameters corresponding to 
Mn atoms in substitutional and T(Ga4) positions. The 
strong deviation of the actual density-functional results 
from this model expectation may indicate that the hy- 
bridization of the Mn d-states with more distant neigh- 
bors has a remarkable influence on the magnetic interac- 
tions. 

The fact that the magnetic behavior of the interstitial 
Mn is rather insensitive to its position in the crystal has 
two important implications concerning the MnGa-Mni„t 
pairs. First of all, the pair as a whole does not interact 
much with the spin of the holes because the effects due 
to MuQa and Muint compensate one another. In this re- 
spect, our ab initio results overcome the troubles of the 



We used the FPLAPW method to obtain total energies 
of supercells simulating various geometric and magnetic 
configurations of Mn atoms in (Ga,Mn)As. In absence 
of other defects, the ground state of the Mn interstitials 
is the tetrahedral T(As4) position. The energy of the 
T(Ga4) position, however, is almost the same. The sit- 
uation changes in the p-type material where the T(Ga4) 
position has a lower energy. The hexagonal interstitial 
position has much higher energy and represents a barrier 
for diffusion of Mn from one to another interstitial site. 
The barrier, and consequently also the mobility of the 
interstitial Mn, depends on the doping. 

The exchange coupling Jpd of Mn interstitials with 
the holes in the valence band is, for both T(As4) and 
T(Ga4), close to the value of Jpd obtained for the sub- 
stitutional Mn. This is not consistent with the simplest 
tight-binding picture of Mn d-states that hybridize only 
with the nearest neighbors. In this way, our result indi- 
cate that the hybridization with more distant neighbors 
may be also important for the magnetic interactions. 

The Mn interstitials are attracted to the substitutional 
Mn and form stable and magnetically inactive pairs. The 
density-functional estimate for both binding energy of 
the pair and for the energy of the antiferromagnetic cou- 
pling is of order of 0.3 eV. This fits well with the present 
day notion of the interstitial Mn in (Ga,Mn)As. [Tll|. 
In contrary to the general opinion, however, we found 
that the efficient pairing is not restricted to Mujnt in the 
T(Ga3Mn) position and we showed the importance of 
the close next-nearest neighbors for the properties of the 
interstitial Mn. 
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